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A3S TRACT 



Optical variability across the con t inentai slope at a 
shelf off Central California was studied using Nimbus-7 
Coastal Zone Color Scanner (CZCS) data. CZCS estimates of 
k(490), the irraiiance attenuation coefficient at 490 nm, 
were expressed as optical depth 1/k(490). A modified atmos- 
pheric correction algorithm was used to account for water 
radiance at 670 nm. Time sequences of 1/Jc (490) were assem- 
bled and partitioned into four zonal transects, at different 
latitudes, spanning Hay through November in 1979, 1980 and 

1982. Empirical Orthogonal Functions (EOFs) were calculated 
for each partition. The first EOFs are dominated by scales 
of order 180 km, with in all cases, a band of low optical 
depth water in the first 100 km adjacent to the coast. 
Scales decrease in successive EOFs, to about 40 km in the 
fifth EOF. The feasibility of joining EOFs from different 
partitions was demonstrated as a precursor for future 
applications to piecewise analysis of oceanic satellite 
da ta . 
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I. 



INTRODUCTION 

Satellite reacts sensing systems offer fast, econcaicai 
means of determining the horizontal structure cf the oceans 
on a global basis. The objective of this thesis is to con- 
tribute to tbs development of empirical methods for using 
satellite images of optical parameters and sea surface temp- 
erature (SST) to infer the upper ocean's vertical structure, 
through interpolation and extrapolation of relatively lim- 
ited in situ data. 

The objective is being approached through regional case 
studies of correlations between optical parameters and phys- 
ical water mass properties in the upper ocean in different 
regions cf the world. Mere specifically, this thesis is a 
preliminary case study cf the California Current region. 
The ultimate goal is to relate statistically the horizontal 
structure cf optical properties observed with the Coastal 
Zone Color Scanner (CZCS) to the underlying vertical struc- 
tures of temperature and salinity, as well as bio-optical 
parameters, for a given region and season. 

The study domain encompasses the continental slope and 
shelf off the coast of California between Point Sur and 
Point Arguellc. This area was selected to investigate an 
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ocean upweliir.g rr 



rhich is known to oersist throughout 



the upwalling season (Traganza,et al., 1979) . The northern 
and southern portions of the study domain are typified by 
complex eddy structure associated with irregular features in 
the bathymetry, such as cff Point Sur. Between Point Sur 
and Point Arguellc, on the other hand, isolines of S ST and 
optical parameters tend to be aligned roughly parallel to 
the underlying isobaths. 

An ensemble of data acquired with the Nimbus 7 CZCS dur- 
ing the summer and fall seasons of 1979, 1980 and 1982 is 

analyzed in this study. The horizontal structure in bio-op- 
tical parameters determined from cloud-free portions of CZCS 
imagery are investigated using a Partitioned Empirical 
Orthogonal Function (PEO?) decomposition. The spatial par- 
titions examined here consist of four zonal transects cross- 
ing the shelf/slope region at different latitudes. The 
specific goals of this analysis are: 

1. To characterize the meridional and- zonal spatial 
correlation structures of ocean color parameters 
(specifically optical depth 1/K(490) in meters). 

2. Tc compare the spatial scales and structures of 
optical variability highlighted by the PEOF 
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decompositions, and to relate these to the historical 
descriptive oceanography of the study region, and 
3. To develop preliminary statistics related to the 
feasibility of joining data from different spatial par- 
titions on the basis of partial subsamples, ana to thus 
provide an optimal interpolation of satellite image 
data into cloud covered areas. 
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II. 



OCEANOGRAPHY CF THE CENTRAL CALIFORNIA 

coisr 

A. THE STUDY DOMAIN 
1 . Cove rag e 

a. Area Domain 

The region investigated in this project is 
located between 32 and 40N, and from the coast of California 
offshore to approximately 126W, Fig. 1. This area was 
selected because it contains water mass structures, includ- 
ing fronts, which strongly influence phytoplankton concen- 
trations, and therefore the optical properties of the ocean 
water. Furthermore, an adequate sample of data was availa- 
ble for this area. 

A subarea of this region is labelled Insert A in 
Figure 1 and presented in greater detail in Figure 2 Insert 
A is bounded by 34 to 38 N, and by 126 to 120w. It is the 
primary study domain of this thesis. The background hydrog- 
raphy and dynamics of this region are described in subsec- 
tions 2, 3, and 4 of this chapter. 
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CALIFORNIA COAST BATHYMETRY 




Figure 1. Ocean Bathymetry Off the California 

Coast (Synt net ic Bathymetric Profiling 
System (SYNBAPS) Data Contoured at 400 a 
Intervals) . 
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b. Time Domain 

In the tine domain, the available CZCS data 
include scenes from summer through early fall seasons in 
1979, 1980 and 1982. Originally, a single season ensemble 
of CZCS data (day through September 1980) was sought. How- 
ever, the time span had to be expanded to three years no 
assemble a reasonably large sample size of cloud-free sub- 
scenes. The sample analyzed contains three scenes in 1979, 
eleven scenes in 1980 and eight scenes in 1982. Detailed 
characteristics of these CZCS images are presented in 
Appendix E. 

2. Geometry 

a. Coastal 

The California coast bounding the area of study 
is characterized by a steep, mountainous coastal range run- 
ning roughly parallel to the coastline. The coastline 
stretching from San Francisco to Point Arguello is oriented 
roughly northwest to southeast, but is interrupted by Monte- 
rey Bay at 36 45'N and by smaller bays in the vicinity of 
Morrc Bay at 35 20' N. No major rivers drain into this 

coastline, although many local rain-generated drainage 
creeks empty hers. 
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Bathyraet r y 



The predominant orientation of the bathymetry is 
roughly northwest to southeast (parallel to the coast) , Fig. 
2. Interruptions of this orientation are evident in the 
vicinity of the Monterey Canyon, Point Sur and the Sur Can- 
yon, the Davidson Seamount, the Taney Seamount and the Santa 
Lucia Banks and Escarpment off Point Arguello. A very 
abrupt shelf break is evident all along this section of the 
California coast. Isobaths tend to diverge south of Monte- 
rey, due to a broadening of the continental shelf and slope 
with distance south, of Monterey. 

3. De scri ptivs Oce an c araph y 
a. Coastal Upvelling 

Coastal upwelling is an oceanic phenomenon which 
has a pronounced impact upon many physical and biological 
processes. Predominantly southward winds during spring and 
summer off the central California coast, yield offshore sur- 
face Ekman transports, which forces compensation water to 
rise from depths of the order of 200 to 300 m (Smith, 1968). 

The upwelling season off the coast of California 
is generally confined to the late spring through early fall. 
The onset of the seasonal upwelling commences in more 
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INSERT A 




Figure 2. 



Ocean Bathymetry Off the California Coast 
(SYNBAPS Data, Contoured at 230 m Intervals) 
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southern waters off the California coast and prograsses 
northward as the season unfolds (Yoshida and Mao, 1957; 
Wooster and Reid, 1963; Pavlova, 1966; Hickey, 1979). The 
CZ CS data set spans the upwelling season and includes images 
from beyond this season into early winter. 

Upwelling has a marked effect on the sea surface 
temperature, causing it to be much lower than would other- 
wise be normal for the latitude and season (Smith, 1968) . 
The relatively lower temperatures are evident in IR images 
of the region (Johnson, 1980; Nestor, 1979). Accompanying 
this decrease in temperature is an increase in surface 
salinity, an upwelling property unique to the regime off the 
west coast of North America (Smith, 1 963) . 

The oceanographic properties of upwelling have 
been documented in many areas of the world, but nowhere with 
the thoroughness of the work off California and Oregon. 
Ship and satellite observations have allowed us to identify 
seasons, centers, and the extent of the upwelling event 
along the west coast of North America. Traganza, et al. 
(1979) used combined satellite and shipboard observations to 

I 

infer nutrient upwelling distributions off the coast of Cal- 
ifornia. Frontal structures and mesoscale eddies that can 
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result from the upwelling phenomenon have been examined with 
relevance to Anti-Submarine Warfare (ASW) by Traganza 
(1979). The use of infrared ( I R ) imagery in the detection 
and description of upwelling was examined by both Johnson 
(1980) and Nestor (1979) . 

The introduction of nutrient rich waters to the 
nearshore euphotic zone greatly enhances the development of 
the in situ phytoplankton population. This enhancement in 
turn causes the upwelled water mass moving offshore to have 
distinctly different optical properties than adjoining off- 
shore waters. The boundary (frontal region) between the 
upwelled water mass and the normal surface water mass is 
thus readily detectable and of great interest. 

Nutrient enrichment off the California coastal 
zone is observed in the regions cf upwelling events. These 
nutrients, which are classified as "biochemically new" on 
the basis of nitrat e-t o-phosphate ratios which approach 
15:1, are brought to the surface from depths up to 300 m. 
By way of contrast, nutrients also present in the open ocean 
surface water approach 5:1 (Nestor, 1979). The added 
nitrates are a primary factor in the increase in phytoplank- 
ton concentrations during the upwelling season. Coastal 
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wafers, on the whole, are fifty times more productive than 
open ocean waters and this difference can be increased dur- 
ing periods of upvelling (Sverdrup, et al . , 1942). The phy- 
toplankton concentrations, with their associated 
chlorophyll-like pigments, have a profound effect on + he 
upwelled radiances measured by the CZCS, as discussed in 
Chapter III. 

Another aspect of up welling, and its relation to 
satellite data, is its effect on regional climate. The rel- 
atively cold sea surface temperature in upweiling zones 
cools the air above and thus increases its relative 
humidity. As a result, low stratus and fog commonly occur 

here in a shallow (marine) layer with warm air aloft. The 
frequent occurrence of low stratus and fog, seriously limits 
infrared and visible satellite coverage during the upweiling 
season. The cool sea water also contributes to a diurnal 
sea breeze by increasing the onshore-offshore pressure gra- 
dient. Onshore winds bring cool, moist air as far as 50 
miles inland (Smith, 1968). 
b. Currents 

The California Current System may be discussed 
and studied in terms of four large scale currents: 
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the California Current, the California Undercurrent# the 
Davidson Inshore Current, and the Sou-hern California Cur- 
rent (Hickey, 1979). The first three of these currents 
directly influence the study domain. Mesoscale currents 
associated with seasonal upwelling are also important hers. 

(1) Cali for nia Current. The California Current 
is a broad wind-driven equatorward current which exhibits 
significant seasonal variations proportional to the changes 
in the wind field (Brown, 1974) . Off Point Conception the 
mean annual location of the current axis is located 270 km 
offshore while the shoreward boundary extends tc 200 km off- 
shore. The current is of the order 700 km wide and flows 
south at 10 to 30 cm per second (Hickey, 1979). 

The California Current is a continuation of 
the West Wind Drift in the North Pacific and flows southward 
along the California coast between 48 and 23N. It turns 
westward between 20 and 30N where it becomes part of the 
North Equatorial Current. This flow regime comprises the 
eastern extent of the anticyclonic NE Pacific Subtropical 

s 

Gyre, which is centered near the Hawaiian Islands (Sverdrup, 
et al., 1942; Chelton and Davis, 1982). 
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(2) California Unde rcurr ent . The California 
Undercurrent, also, referred tc as the California Counter- 
current, is the poleward subsurface flow over the continen- 
tal slope. Maximum poleward flew occurs during the summer 
and fall seasons at depths of 200 to 250 m (Pavlova, 1966 
and Hickey, 1979). The flew can be described as a broad 
current with a central jet. It is this jet structure that 
is most often measured and referred to when applying 
specific values to the Undercurrent. The broad poleward 
flow has a gecstrophic component alongshore near the shelf 
break of approximately 15 cm/sec (Ccddington, 1979). 

The flow appears to have a jet-like struc- 
ture, both vertically and horizontally, and to extend to the 
bottom over the slope. The existence of a high speed jet 
core of the order of 20 to 70 km in width, was first sug- 
gested by Beid (1962, 1963). Subsequent direct measurements 
of these jets have produced values as high as 40 cm/sec off 
Northern Eaja and values of 16 cm/sec off Washington 
(Wooster and Jones, 1 970) . The depth of the high-speed cere 
varies seasonally. It rises from depths of 200 to 300 m to 
the surface during the late fall and winter north of Point 
Conception. Here in is referred as the Davidson Inshore 
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Current by 31 any authors (Hickey, 1979; Pavlova, 1966; Ingra- 
ham, 1967). Event-scale fluctuations (of rhe order of ICO 
km and 10 days) in the flew appear to be correlated with the 
alongshore component of wind stress (Nelson, 1977) . 

The extent and time scale of continuous 
alongshore flow, and the width of the region of northward 
flow below 500 m, are important topics yet to be answered 
about the California Undercurrent (Hickey, 1979). 

(3) iDE^cre Current . North of Point 

Conception, the poleward surface flow in the nearshore 
regions off the West Coast is known as the Davidson Inshore 
Current. It is associated with winter weather circulation 

* 1 

patterns. As the southward winds weaken and tend toward a 
northwestward flow, the Davidson Inshore Current becomes 
established (Hickey, 1979) . The current flows near the 
coast, usually within 100 km, well inshore of the California 
Current and is confined to the continental shelf and slope. 
Pavlova (1966) reported that north of Point Conception, the 
Davidson Inshore Current reaches its maximum development at 
depth (200 to 250 m) in the summer and autumn. In August 
the Davidson Inshore Current is scarcely noticeable at 
the surface despite active development a depth. Maximum 



27 



surface development is reached from October through April, 
i. e. , late autumn to early spring. In December, the core of 
maximum velocity emerges at the surface ana nr. the late 
spring it almost completely disappears (Reid, 1960; Reid, 
et al., 1958; Pavlova, 1 S66) . Poleward velocities of up to 
25 cm/sec were recorded (Reid and Swartzlose, 1962) within 
80 km of central California in January. 

The Davidson Inshore Current and the Cali- 
fornia Undercurrent are often discussed as though they were 
separate currents. Both currents transport Equatorial- type 
water northward at least as far as Cape Mendicino (Pavlova, 

1966). Also, no subsurface maximum has been found in the 

• * 
flow of the Davidson Inshore Current. These characteristics 

support a view that the Davidson Current is simply the sur- 
face expression of the California Undercurrent, rather than 
a separate current superimposed on it. 

(4) Other Currents. The presence of eddies 
throughout the the California Current System has been docu- 
mented for many years (Bernstein, et al. , 1977). The time 

scales for these eddies, as well as the processes responsi- 
ble for their generation and subsequent dissipation, is an 
area of increasing study. 
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Between San Francisco and a point about 
half the distance to Eoinz Conception, there is a permanent 
counterclockwise eddy than produces northward flow during 
all months except April (Brown, 1974). A second eddy just 
north of Point Conception forms during the summer months and 
makes northward flow continuous from Point Conception to San 
Francisco (Hickey, 1 979) . Willmott (1983) has shown that 
these features may be produced by flow separation of the 
California current in the vicinity of major coastal capes. 

Reid, et al. (1963) made direct measurements of 
an eddy (90 km in diameter) off the northern coast cf Baja 
California. Hypotheses for eddy formation discussed in 

I 

their paper are as follows: 

(1) The process of upwelling and the offshore movement of 
the colder, more saline waters might degenerate into eddies. 
The lateral shear between the upwelling flow away from the 
coast and the California Current and Undercurrent (bare- 
tropic instability) could produce eddy structures. Tempera- 
ture and salinity differences set up strong baroclinic zones 
/ 

along the upwelling boundaries which could result in eddy 
formation. (Sverdrup and Fleming, 1941) 
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(2) The offshore surface flew during upwelling should pro- 



duce a counter current (Munk, 1950) . If there is 
substantial nerrh-south variation in the intensity of the 
winds, then seperate countercurrents of different strengths 
might occur along the coast. 

(3) A second hypothesis proposed by Reid, et al. (1963), is 
that the deeper countercurrent may transfer momentum upward 
to the surface layers, at times when, or in regions where, a 
surface current does not prevail. This could cause spot 
intrusions of colder circulating waters that form eddies 
where neither surface countercurrents nor coastal upwelling 
produce them. 

i 

Additionally, the effects of bottom steering by coastal 
topography, and the associated trapped motions must be con- 
sidered when discussing eddy formation . (Hurlburt, 197U; 
Johnson, 1982; Willmott, 1983) 
c. Water Masses 

Descriptions of the water masses that contribute 
to the California Current System are given in Tibby (1941), 



Sverdrup, et al., (1942) and Reid, et al., (1958). Four 
major sources are discussed by the authors: 

(1) Subarctic Water Mass - from the north. 

(2) Central Water Mass - from the west and northwest. 
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(3) Equatorial Water Mass - from the south. 

(4) Water derived from upweiling sources. 

These sources were simplified in Tibby (1941) and Sverdrup, 
et al., (1942) into two extreme sources named "Subarctic 
North Pacific" and "Equatorial Pacific". 

The percent of each water mass comprising a sam- 
ple can be defined by entering Figure 3 with a T-S pair. 
However, the determination of percentage composition by this 
means cannot be used for water above depths of about 100 m. 
This restriction is due to vertical mixing in the nearsur- 
face layer related to the effects of wind and local changes 
due to hear and mass fluxes across the air-sea interface. 
Any mixing along surfaces of constant a t would be severly 
masked in these shallower depths by the effaces of turbulent 
vertical mixing. Also, below 1000 m the differences in the 
T-S relationships of the two extreme water masses are negli- 
gible. For intermediate depths, as might be expected, the 
percentage of equatorial water decreases in the direction of 
northward flow. The Undercurrent is characteristically 
warmer and more saline than the California Current, and it 
has a salinity maximum cn the a = 26.54 surface. Cff 
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Monterey and below 3 00 m, the water is greater than 6 3 1 
Equatorial Water and this percentage increases both with 
depth and movement towards lower latitudes (Brown, 1974). 



i 
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Figure 3. 



Graph Showing T-S Curves Defining 
and Equatorial Pacific Water, ana 
Various Percentages of Equatorial 
Assuming Mixing along Surfaces of 
(Brown, 1974). 



Subarctic Water 
Curves for 
Pacific Water 
Equal a 
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A. INTRO E OCT ION 

The physical processes of absorption and scattering 
relate the upwelling radiance just beneath the sea surface 
to the constituents of the water (Gordon, 1976) . Except for 
coastal waters and waters influenced by river discharge, 
biological constituents play a dominant role in these pro- 
cesses (Smith and Baker, 1978; Jerlov, 1976). Optically, the 
most important biological constituent is phytoplankton, 
microscopic plant organisms that photosynt he size and make up 
the first link of the oceanic food web (Steele, 1970). 
Ch lorophy 11- a‘ is the dominant p hotosy nthe tic pigment, and 
absorbs light strongly in the blue and red regions of the 
visible spectrum (400 to 700 nm) (Hovis, et al., 1980). 
Therefore, as the concentration of phytoplankton increases, 
the color of the water is shifted toward green hues from the 
deep blue of its pure state. By measuring upwelled radiance 
(backscatt ered daylight) in specific spectral bands, we can 
determine the concentrations of phytoplankton pigments in 
the ocean (Gordon, et al. , 1980; Gordon, et al. , 1983). 
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This chapter first describes the CZCS sensor and ins 
capabilities, ar.d then the measured signal is discussed. 
Algorithms that are currnstly applied to this signal to cor- 
rect for atmospheric effects are discussed. Finally, the 
algorithms designed to convert the corrected radiance values 
to phytoplankton concentrations, C, and irradiance attenua- 
tion coefficient, k, are presented. 

B. SYSTEM DESCRIPTION 

1. The N imbus- 7 Coas t al Zone Colo r Sc anner (CZCS) 

The CZCS was built by the Ball Brothers Research 
Corporation to NASA's specifications. The instrument is a 
spatially imaging multispectral scanner. Six spectral bands 
are precisely coregistered and internally calibrated. The 
swath width of the CZCS is slightly more than 1600 km. 
Characteristics of its five visible (443, 520, 550, 670, 750 
nm) and one thermal IR (10.5 to 12.5 ym) channels are summa- 
rized in Table I. The CZCS has an active scan of 78 degrees 
centered on nadir and a field of view of 0.0485 degrees, 
yielding a geometric instantaneous field of view of 825 m 
(at nadir) from a spacecraft altitude of 955 km. It can 
tilt the scan plane 20 degrees from nadir in 2 degree incre- 
ments along the satellite track to minimize the influence of 
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direct sun glint. The Nimbus- 7 spacecraft is in a sun- 
synchronous orbit with ascending nods near local neon. 



TABLE I 

Char actsr istics of the CZCS 
(Hovis, et a 1. , 1980) 



3and 

Number Wavelength Gain 

(nm) 

3 

1 433 to 453 2 

1 

0 

3 

2 5 10 to 530 2 

1 

0 

3 

3 540 to 560 2 

1 

0 



4 660 to 680 2 

1 

0 

5 700 to 800 

6 10,500 to 12, 500 



Saturation 


Measured 


Radiance 
(mW/cm 2 sr ym) 


signal/noise 


5.41 


7.64 

9.23 

11.46 


158/1 


3.50 


5.10 

6.20 

7.64 


200/1 


2.86 


4.14 
5.10 
6 .21 


176/1 


1.34 


1 .91 
2.32 
2.88 


118-/1 


23.90 


350/1 
0.22 K* 



* Noise equivalent temperature difference at 270 K. 
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The designed purpose of the CZCS experiment was to 
provide estimates of the nearsurface concentrations of phy- 
toplankton pigments (defined to be chlorophyll-a and its 
associated degradation products, called 11 phaeop igment s ") by 
measuring the spectral radiance backseat tared out of the 
ocean (Gordon and Clark, 1581) . The radiance scattered out 
of the ocean that reaches the sensor is a very small portion 
of the total radiance received. Consider the physical set- 
ting where solar irradiance F„(X) at a wavelength X is inci- 
dent on the top of the atmosphere at a zenith angle 0 Q and 
azimuth 0 q and the scanner is detecting total radiance L C (X) 
at a nadir angle 0 and azimuth angle 0 . L t ( X) consists of 
radiance which has been scattered by the atmosphere and sea 
surface, radiance generated by Fresnel reflection of the 
direct (unscattered) solar irradiance from the rough ocean 
surface (sun glint), and solar irradiance scattered from 

beneath the sea surface t,(X)L (X), where t ( X) is the diffuse 

d w d 

tr ansmitt ance. 

Observations (Gordon, et al., 1983; Gordon, et al. , 
1980) produce values cf L t (X) in the blue that are ten times 
greater than L (X) . These effects are principally due to 
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scattering by the air (Rayleigh scattering) and by micro- 
scopic particles suspended in the air (aerosol scattering) , 
both cf which increase the radiance detected at the sensor. 
Fresnel reflection (sun glint) can be ignored as the tilting 
capability of the CZCS minimizes its effect. However the 
scattering effects, both Rayleigh and aerosol, must be 
removed from 1 1 (X) to give usable values for the upwslled 
radiance L (X). 

W 

C. CZCS GEOPHYSICAL ALGORITHMS 
1. At mo s ph e ric Corrections 

From the signal description in section A of this 
chapter, we can construct the following formula 



L t (X) = L r (X) + L a (X) + t d (X)L w (X) 

where 



( 1 ) 



L t = Total radiance 

L r = Radiance due to Rayleigh scattering 
L a = Radiance due to aerosol scattering 
L w = Upwelled radiance from beneath the sea surface 
t d = Diffuse transmittance of the atmosphere 
X = fcavelength 
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A S 



the total 



previously mentioned L t is 



a diancs 



measured by the CZCS. The Rayleigh scattering term can be 



expressed as 



L r (A) 



F q (A) x r (A) 




+ { p ( 0 ) + p(0 Q ) >P r (a + ) 




COS 0 



Where 



F = The instant aneo us extraterrestial solar 



o 



irradiance. 

t r = The Rayleigh optical thickness of the atmosphere. 

F = The Rayleigh scattering phase function. 

a_ = The scattering angle through which photons are 
backscatter ed frcm the atmosphere to the sensor 



without interacting with the sea surface. 

a+ = The forward scattering angle of those photons 
which are scattered in the atmosphere toward the sea 
surface (sky radiance) and then specularly reflected 
frcm the surface into the field of view of the sen- 
sor ( p ( 0 ) term) as well as photons which are first 
0 

specularly reflected from the sea surface and then 
scattered by the atmosphere into the field of view 
of the sensor ( p ( 0 ) term). 

P = The Fresnel reflectance of the air-sea interface. 
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T = The twc-way ozone transmittance of the ataoschere. 
0 3 

0 = The sensor zentith angle at the observed point on 

the sea surface. 

9 = The solar zenith angle at the observed point on 

o 

the sea surface. 

The aerosol scattering term is found using the 670 
nm channel, where there is only a negligibly small contribu- 
tion by the L w term. (This is referred to as the "black 
ocean" assumption.) He calculate aerosol radiance at X = 
670 as 

L ( 670) = L ( 670) - L (670) . (3) 

a. C L 

The key assumption in this algorithm is that the ratio of 
aerosol wavelengths is constant over a scene, and is given 
as 

L (X) 

£ ( X , 6 7 0 ) = — ‘ (4) 

L a (670) 

e is calculated using either simultaneous direct radiance 
measurement from a ship, or upwelled radiance values mod- 
elled at a clear water pixel (Gordon and Clark, 1981) . The 
latter method is discussed in section 2 of this chapter. 
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5 



Returning to equation { 1) 



d 



(A) is the diffuse 



transmittance of the at. me sphere and sea surface, which ray 



be approximated as 



[1 - P ( 0 ) ] [ (x r (X) 



+ x „ (An/ cos 6 
0 3 



(5) 



t d ( A) 



exp - 



where all tens have teen previously defined except m , which 
is the index of refraction of water relative to air and is 
assumed to be 4/3 for the wavelengths {400 - 700 nm) . 

We have now developed the basics for extracting the 
upwelled radiance values, from the CZCS detected signal. 



2. Clear Water Rad ia nee 

The scene constant, e , given in equation (4) is cal- 
culated as 



L 



t 



£ ( A , A ) = (A/A ) n 

0 0 




( 6 ) 



F (A ) exp [ — x _ (A )(sec0 + see© )] 
oo 0 3 o o 



where n is called the Angstrom coefficient 



Equation (6) 



can be rewritten as 



s ( A , A ) = (A/A ) n 

O O 



F (A) T (A) 
o 0 3 




(7) 



where T Q ^ ( A) = exp|^-x 0 (A ) (sec 0 + secB, )J and where for 
the CZCS A q = 67 0 nm . 

Gordon and Clark (1981) developed the concept of 

clear water radiance for atmospheric correction of CZCS 

imagery. The strategy employed in this study was to find an 

area cf the image that cculd be assumed to have a chloro- 

3 

phyll concentration less than 0.25 mg/m . At this low con- 
centration, L at 520 and 550 nrn are assumed to be 

W 

essentially constant for a given solar elevation. Then, 
given these “clear water values” cf L W (A) at one position, 

i 

L (A) as calculated using equation (1). e(A,670) are found 

from equations (4) using the computed L ( A) and L a (670) 
value. Finally, rearranging equation (7) we find that 



n ( A) = 



) e ( A , A ) / 


f F (A) T „ (A) 
o 0 3 


j 


1 0 


F (A ) T (A ) 

L o o O 3 o J 


j 



( 8 ) 



In ( A / A ) 
o 



Values for n at 520 and 550 nm ( n (520) and n(550)) are com- 
puted, then averaged to estimate n(443). The Angstrom coef- 
ficient at 443 nm cannot be directly measured in this way 
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because L (442) is highly sensitive to even minute fluctua- 

W 

rions in chlorophyll concentrations at low concentration. 

An important aspect of this algorithm is that nei- 
ther surface measurements of L w (A) , nor any properties of 
the aerosol are required to implement the atmospheric 
correction . 

3. Pa ram eters 

a. Chlorophyl Concentrations 

Determination cf chlorophyll concentrations 

C frcm ratios of L (A) relates the surface value of C to the 

w 

ratio of the upvelled radiance at two different wavelengths 
(Morel and Prieur, 1977; Gordon and Clark, 1981). The basis 
for this is that to a first approximation L is proocrtional 

W . 

to the ratio cf the volume backseat tering coefficient, 

8(A) b (A) , and the volume absorption coefficient, a (A), of 
the water plus its constituents (Gordon, s: al. , 1983) . The 
contributions from the individual constituents can be summed 
to provide a total value for each optical coefficient. 
Moreover the contributions to B(A), b (A) , and a (A) arising 
from phytoplankton and their pigments are assumed to be pro- 
portional tc chlorophyll concentration C. Taking a ratio of 
L w at two different wavelengths and applying the assumption 

L (A) a B (A ) b (A ) /a (A) , we obtain 
w 
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( 3 ) 



L ( A i > B ( A i ) b ( A i ) a ( A 2 ) 

w = j ; 

L^ ( A 2 ) B(Aj)^Oo^ 

3= cause of -ha non-linearities involved in -he individual 

constituent contributions to B (A ) , b (A ) , and a (A ) , we appeal 

to a heuristic observation that 



L 

w 



( A i ) 



L 

w 



<A 2 ) 



R(C,K, . . . ) 



( 10 ) 



i. e. , the ratio of two upwelled radiances is a function R 
of the chlorophyll concentration, C, and the diffuse attenu- 
ation coefficient, k, as well as ether optically important 
constituents of seawater. It was then assumed that 3 is 

t 

related to C through a log linear model of the form 



Log C = Log A q + A i Log R ( A i , A 2 ) » (11) 

which was empirically fit to observations to determine coef- 
ficients A, and Ai . Thus, pigment concentrations are com- 
puted from CZCS data using the sguation 

A i 

C = A R (12) 

o 

The empirical coefficients presented by Gordon, 
at al. (1983) have been adopted by NASA and are: 

Case I: C < 1.5 for R (443,550) 

A = 1. 1297959 

o 

A ! = -1.71 
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Case II: C > 1.5 for R(443,550), bat 

C < 1. 5 for E (520, 5 50} 

A = 1.1297950 

o 

Ai = -1.71 

Case ill: C > 1.5 for E (5 20, 550) 

A = 3. 3265955 



Ai = -2.44 

3 

where C is ip mg/m . 

b. Diffuse Attenuation Coefficient 

A similar development of the algorithm for the 
determination of the diffuse attenuation coefficient, k (A ) , 
is given by Austin (1981) . Like the chlorophyll concentra- 
tion algorithm, this algorithm derives a value based on the 

ratic of L at two wavelenaths. k can be defined as 
w 



k(A) = 



-1 



dF( A , z) 



F ( A , z ) dz 



(13) 



Equation (13) can be solved for irradiance F (a ,z ) to obtain 



F ( A , z 2 ) = F ( A , z ! ) exp £-k(A)(z 2 -Zi)J 



(14) 



He nee 



k ( A ) = 



z 2 -z l 



In 



F ( A , z 2 ) 
F ( A , z ! ) 



(15) 
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Empirically derived coefficients 



from spectral data yield 



k(490) = 0.0833 

(Austin, 1981). 

D. SIGNAL FACTORS 



L (443) 
w 

L (550) 
w 



1 . '4 9 1 



+ 0 . 022 



m 



( 16 ) 



Many factors have been accounted for with these algor- 
ithms by either mathematical and empirical models or heuris- 
tic assumptions. The determination of the total radiance 
values in the first four channels of the CZCS allows us to 
apply the corrections to determine upwelled radiance. The 
constituents of the water which affect its absorption and 
scattering properties are then empirically derived. 

The distribution of phytoplankton is controlled by many 
local, mesoscale and global factors, including solar radia- 
tion, global weather patterns, and ocean circulation pat- 
terns. The mesoscale events of upwellirg or eddy 
circulation can have important regional effects. These fac- 
tors are too numerous and varied to be modelled on a theo- 
retical basis. However, empirical modelling can produce 
relatively accurate and consistent results. 

The measurement of these bio-optical parameters from 
space allows us to remotely determine their relationships to 
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physical events in t he regions under study. Time scales, 
spatial scales and specific features can be discerned using 
-he known (or hypothesized) relationships between inherent 
optical properties of the ocean water constituents and the 
forcing involved in their distribution. 
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THG3QNAL FUNCTION ANALYSIS 

IIH2225 



A. INTRODUCTION 

The concept of principal component analysis has been 
presented and utilized in different forms over the past 
eighty years. Fitting a line tc a data set was usually 
accomplished using a least squares method. Distances to 
this line from each point were measured parallel to an arbi- 
trarily set axis. From the early work of P earson ( 1 90 1 ) , 
this method was adapted sc that the perpendicular distances 
from each point to the best fir line were measured. Figure 
4 illustrates this difference and shows that the first 
method is tied to a coordinate system while the Pearson 
approach is independent of coordinate systems. This new 
method laid the foundation for the development of principal 
component deccmposi ti on techniques. These techniques have 
since been utilized in many forms and referred to by similar 
names in a number of disciplines. Applications in psychol- 
ogy by Eckert and Young ( 1936,1 939) , although somewhat dif- 
ferent in their development, contain the essential elements 
of data analysis and principal component decomposition as 
used in geophysical disciplines today. 
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Figure 4. 



Plot Showing the Difference Bet we 
of Distances to a Line Parallel t 
Values) and Perpendicular to the 
(Preisendorf er , et al . , 1980) 



sn Minimization 
o an Axis (d 
Line ( p Values) . 
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Loren z ( 1 955) , 



Meteorological applications by 



Vn 4- r- — 



bach (1967) , and Rinne, et al. (1979) demonstrated the convie- 
nence of representing a large climatological field with a 
smaller set of values. These areas include: 

(1) Non-linear statistical prediction (Lorenz, 1 956 ), 

(2) Ncn-linear dynamical prediction (Lorer.z , 1956) , 

(3) 500 mb height field representation (Rinne, et al. r 
1979) , and 



(4) Sea level pressure, surface temperature, and pre- 
cipitation pattern representations (Kutzbach, 1967) 
Other uses of EOF analysis techniques in oceanography 
include the representation of ocean color spectra (Mueller, 
1976) and of wave spectra (Aranu vachu pun and Thorton, 1983). 

The principal difficulties encountered in principal com- 
ponent analysis problems relates to the selection of the 
'meaningful' subset of components and to oheir physical 
interpretation. Methods of selection of the principal com- 
ponents are also widely varied. Preisendorf er , et al. (1981) 
discussed two methods which together involve seventeen dif- 
ferent testing rules. Empirical selection of a cutoff value 
for variance or forcing factors can also be utilized. 



50 



Visual inspection of the data which leads to a clear cut 
(albedo subjective) choice is alsc an option. 

The enormous dare volume inherent in satellite data sets 
begs application of the techniques of principal component 
analysis. Principal component analysis techniques often 
allow the efficient representation of a large data set by 
its first few principal components with a negligible loss of 
information. The advantage gained is reduction in the num- 
ber of variables needed tc represent the data. Reducing a 
data set to its principal components can also aid in the 
interpretation cf the data by separating noise from the sig- 
nal. Principal component analysis theory can be applied to 
preliminary explorations within a relatively unstructured 
domain of knowledge, one in which the fundamental laws gov- 
erning the processes under study are still being 
defined. (Preisendcr f er, et al., 1981) 

A brief review of the EOF analysis follows to provide 
background for the later analyses. The reader is referred 
to Pries endor fer , et al. (1981) for a more complete devel- 
opment and history. The following matrix algebra notation 



is adopted throughout this thesis. 

1. No underscore denotes a scalar X 

2. A straight line underscore denotes a vector .... X 
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3 



. A curved lire underscore denotes a matrix ..... 

4. A straight line cverbar denotes a mean value... X 

5. The use of a superscript "T" denotes a matrix 
transpose. 



8. EOF EQUATIONS 

1. Raw Data Con v er sic n 

Following Preisen dorf er , ex al. (1981) let be the 
raw (uncentered) data matrix. 



f ' (1,1) f ' (1 ,2) ... f ' (1 ,p) 

f ' (2,1) f ' (2 ,2) ... f ’ (2 ,p) 



F ' 



f ' ( n , 1 ) f ' ( n , 2 ) 



f ' (n , p) 



(17) 



where f (i, j) is the measurement in the i ' th time point and 
j'th spatial position. In the present investigation, each 
member of F’ will correspond to an optical parameter meas- 
ured by the CZCS at a particular time and spatial position. 
To convert the raw data matrix, F', to a centered data 
matrix, F, the temporal means are computed and subtracted 
from JF' . The temporal mean vector f(x) is calculated as 



f (x) 



i n 

1 T. 

= n 



f ' (t.x) 



( 18 ) 



t = 1 

The centered data matrix, F, is then defined as 
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( 19 ) 



— : 

f * (1,1) - f(l) f ' (1,2) - f (2) . . . f ' (1 , P ) - f(p) 
f 1 (2,1) - f ( 1 ) f ' (2,2) - f (2) . . . f ' (2,p) - f(p) 



F = 



f ’ ( n , 1 ) - f(l) f ' ( n , 2 ) - f (2) . . .f ' (n,p) - f(p)_ . 

Each element f(t,x) of J consists of a raw data measurement 
with the temporal mean removed. The centered data matrix , 
F, can be written as 



F = F ' - F 



( 20 ) 



where £ is the matrix containing as rows the transpose of 
the mean vector f T = £ f ( 1 ) , f (2) , . . . , f (p) J . 

2. Pri ncipal Direction of Scatter 

To find the direction, ej, (in the space domain) 
along which the scatter (or variance) of the data set is 
greatest, consider the projection of the data vectors f_(n) 
along an arbitrary direction e^ 



D( t , e ! ) = f T ( t ) e l 



( 21 ) 
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Squaring this length and summing over all n observations, 
aives a measure cf the scarier cf the data along the direc- 



tion, e i , namely 



2 n 

D (ej) = Z 

t = 1 



T 

f (t) 






(22a) 



The righthana side cf (22a) can be expanded to yield 
n 



2 

D (e^ = 



-2 ^cofco.J 

-sITe 5l • 

The next step is to define the "Scatter Matrix", S, 



(22b) 

(22c) 



S = F F 

/v 



(23) 



with elements 



Expanding 



n 



s(i,j) =2 f ( t ) f ( t ) 

t = r 

the above equation produces 



(24) 

for each member of S 



n 

s ( i , j ) = 2 ( f ' ( t , i) -f ( i) ) ( f ' ( t , j ) -f ( j ) ) T . (25) 

t = 1 

If the matrix is normalized by dividing by (p - 1) , then 

when i/ j, the members of are covariance values, and when 
i = j (the trace of the matrix) the members are variance 
values (i.e. each element is the variance of f at a single 
spatial grid point) . The scatter matrix, S, is symmetric. 
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Therefors, it generally has p no n-negative eigenvalues 1 (j 
= 1 f p) and associated eigenvectors e. (j = 1, p) , 

provided that the rank of is equal to p. 

3* Prin cip al Component, Big enval ue and Eige nvec tor 
Rep rese nt a t en ” ” ” “ 

The first principal component of an observation vec- 
tor f is defined to be the linear combination 



a l = e l l f l + e 1 2 ^ 2 + •••+ e i p f p = £ I 

whose sample variance 
P P 



( 26 ) 



; 1 S e i 

1 V'l - 1 



( 27 ) 



S a ' Z) S ei i e lj S ij ' 

1 i = 1 j = 1 

is a maximum for all possible vectors e, subject to the con- 



straint that 



T i 

eiei = 1 



( 28 ) 



Introducing the Lagrange multiplier X i, the maximum variance 
must satisfy 






s +X i ( 1-e i e i 
a i 1 - L - L 



"177 + A 1 < i-eTe 1 )J 



= 2 (S - X i I ) e i = 0 

t/i ^ - 



( 29 ) 



For nen-trivial solutions , X lr must be chosen such that 



5 " 



= 0 



( 30 ) 
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and, therefore, Xi , is an eigenvalue of S, and 3 [ is its 

u' *** 

associated eigenvector. Furthermore 



b e 1 

'V. - 1 



= X 



and since e ie i = 1 



T 

e S e = X 
“i^ - i l 



2 

s 

a 1 



(31) 



(32) 



i.e., the first eigenvalue of ^ is int erpr etable as the sam- 
ple variance of S. If we expand this development to the 
other eigenvalues and eigenvectors of S, we obtain 



E = j'e.j ; j = 1 , . . . p , 
the eigenvector matrix and 




the diagonal eigenvalue matrix. 



(33) 



(34) 



Thus , 

s \ - i \ • < 35 ' 

T T 

In terms of E the constraint e. e = 1 becomes E E = I. Where 
I is the identity matrix. Therefore, if we multiply both 

T 

sides of eguatron (35) by E * we obtain 



S 



T 

E L E 



'x, 'x, 'x. 



(36) 
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Now, using the definition provided for f and equation (20) , 
the principal component matrix can be defined as 



$ - i i . < 37 > 

This is the desired principal decomposition of ? where 



= (F 


E) 


T 

1 l 


(38a) 


T 


T 






= E 


F 

Ow 


F E 


(38b) 


= S T 


S 

% 


E 


(38c) 


- f 


E 

% 


L 


(33d) 


= L 




. 


(3 8e ) 



C. PARTITIONED EOF ANALYSIS 
1 . P urgose 

The EOF analysis method outlined above works very 
well for a large continuous data set. However, geophysical 
data sets are rarely continuous. In the case of satellite 
data, cloud cover results in many gaps. Sometimes these 
gaps can he bridged by linear interpolation, e.g., when they 
are small and surrounded by good data. Often this is not 
the case and so a scheme of utilizing non-con tinuous data is 
necessary. Here, the purposes of partitioning are: 
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(1) To maximize the sample size in she presence of 
cloud cover, thus allowing statistical computations for 
subregions ; 

(2) To highlight spatial structures of variance fea- 
tures locally, before absorbing them into the modes of 
the overall domain; and 

(3) To achieve computational conviencs. 

Briefly, partitioning permits EOF analysis using 
small subsets of the overall data set. These subsets are 
partitioned to yield continuous data in each subdomain. An 
EOF analysis is completed on each individual subset, and an 
eigenvalue matrix, an eigenvector matrix and a principal 
component matrix are obtained. The next step is to perform 
an EOF analysis to join the principal components of the sub- 
sets. This second EOF analysis produces 'joining functions' 
which relate two non-overlapping subsets. 

2. Rules and Methods 

When performing the partition of any data set cer- 
tain rules must be observed to maintain the statistical 
reliability of the computations. Two obvious and basic 
rules are; 



58 



(1) The minimum partition size (number of pixels} aust 
be greater than or equal to the sample size (i.e., if 
there are 25 sample days each partition must have 25 or 
more pixels). In practice, the spatial dimension will 
he required to be significantly greater than the sample 
size . 

(2) The partition size should net be so small that the 
spatial structure is dominated by noise (e.g., a parti- 
tion boundary will not be placed in a major feature, 
such as a front or eddy of length scale much less than 
the partition size). 

The methods involved in the partitioning are subject 
to the abeve principles, together with a general understand- 
ing cf the physical processes occurring in the study domain. 
Four tracklines at 35N, 35 22'N, 35 40' N and 35 53' N, were 

used to aid in this initial trial of partitioning (Fig. 5). 
The radiance values and computed optical parameters along 
each trackline were plotted versus distance from the coast. 
These plots were aligned to pictcrialiy represent, the data 
and its gaps (due to clouds). The partitioning scheme was 
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Figure 5. 



Trackline plots. 
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then applied to try and produce subsets that were as com- 
plete (continuous) as possible over the time domain. The 
total data matrix is thus partitioned into P subdomains 

i t - i p ■ >> - • 1391 

where the subscript t denotes the total data set and sub- 
script p denotes the partitions of the data set. Each F is 
the data matrix for grid points falling within grid parti- 
tion p, and contains all time points for which complete data 
were acquired in that subdomain. 

3. Eg ua t ion Development 

The EOF decomposition discussed in B, is applied to 
each partition separately, such that for each partition, p, 
the scatter matrix is given by 



and from (36) 




(40) 



S = E L E T (41) 

a, p p p 'v, p , 

where jS = Eigenvector matrix of the spatial partition p, 
and 



L 

oO P 



Eigenvalue matrix for partition p. 
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The matrix of principal components for each partition is 
given by 



A = F E 

^ p ^ p 



Equations (38) require that 



L = A A 

'Vr p *\J p ^p 



and so (41) can be written as 



T T 

S = E A A E 

^P ^p ^P 



(42) 



(43) 



(44) 



Relating this to a global scatter matrix, 

(45) 

Where <^i 2 represents the matrix of covariances between grid 
points in domains 1 and 2, and so forth. 

Now the joining process is developed. For any num- 
ber cf partitions (two are used in this development). 
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where the subscripts denote partition number. 

For this combined set the scatter matrix is giver, by 
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Using equations 
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The joining functions* J, are defined as the eigenvectors of 
the central matrix given in (49) . Finally using (4 1) , (43) 

and (44) , (49) becomes 
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(50) 



The joining functions, J, relate the separate subscenes to 
each ether across an overall study domain. The interpreta- 
tion of these functions should allow examination of the var- 
iations that occur throughout the domain and localized 
affects on the individual partitions of the domain. 
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D. 1ST EE PE STATION 

The coordinate system defined by the eigenvectors gives 
the domain for the principal components. In the present 
study, as well as most geophysical applications, the princi- 
pal components can be thought of as temporal amplitudes and 
the eigenvectors as their spacial modulators (Preisendorf er, 
et al., 1981) . 

The i-th principal component is that linear combination 
of the data field which explains the i-th largest portion of 
the total field variance. Essentially the eigenvectors 
define a direction of variance, while the principal compo- 
nents give the amplitude of the variance in the direction of 
the associated eigenvector. 

Once a data sat has teen reduced to a set of eigenvec- 
tors and associated amplitudes, the question of signal ver- 
sus noise arises. A decision as to which components of the 
data field have significance, and which components of the 
data field have no physical meaning must be made. Some sort 
of a selection process must be defined and applied. Bases 
for these selection processes should have their roots in the 
physical processes being studied. 
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A. INTRODUCTION 



RESULTS 



The major focus of This thesis was to achieve a first 
step towards the analysis cf the obtained data set. Much of 
rhe energy in producing these results was directed toward 
the processing of the data to a usable form for EOF analy- 
sis. Appendices A and B give a detailed discussion of the 
processing techniques utilized and an accounting of all 
adjustments applied to the data. 

The results obtained in this thesis encompass three dis- 
tinct areas. The first result emerged from the data 
processing and a discovery of the breakdown in the black 
ocean assumption. The remaining areas are interrelated as 
one deals with the data set prior to the EOF analysis, while 
the ether attempts to relate this to a statistical meaning 
using EOF analysis methods. 

B. CORRECTIONS FOR NON-ZERO L ( 670) IN COASTAL WATERS 

w 

Preliminary examination of this data set showed that, 
near the California coast, the assumption that L w (670) = 0 

breaks down (Chapter III, Section C) . This finding pre- 
sented a need for an adjustment algorithm. 
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At pixels where the upwslled radiance at 670 r.m , 
L (670), is significantly greater than zero, calculated val- 

W 

2 

ues cf L (443) ars oftan < 0.01 (raw/ (era " -sr-m) ) (approxi- 
w 

mately 1 digital count in CZCS channel 1). This is 
unreasonable even in moderately turbid ocean waters. Smith 
and Wilson (1981) observed that in coastal waters off Cali- 
fornia, where pigments and/or sediment concentration are 
relatively high, it is not uncommon for the subsurface 
upwelling radiance L w (670) to be non-zero. They developed 
an iterative procedure to account for this, which is similar 
to that developed independently and used in the present 
processing . 

The procedure involves two major steps. 

In the first step, which is invoked when L w (443) < 0.01: 

1. Set I (443) = 0.01 , a minimal value for daylight 

W 

tackscatter and slightly less than one digital count in 
CZCS channel 1 . 

2. Decrease L (670) and increase L„(670) to be consis- 
tent with the new value of L (443) (using equations (4) 

W 

and (1)). 

3. Recalculate L (52C) and L (550). 

w w 

4. Recalculate C i, C 2 , K (4 90) , K(520). 
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The second step is based on the assumption that the 
C 2 algorithm (equation (12), case III) is robust and insen- 
sitive to moderate errors in L (670) . This assumption was 

a 

supported by sensitivity calculations which showed C 2 values 
to vary by less than 30 percent for wide variations in e (X 
,610). If C 1 - C 2 > 0-5 and C 1 > 2, the correction is pre- 
sumed unreasonable due to L (670) being too large (and cor- 
respcndingly , L (670) being too small). Next, values of 

w 

K(490) and K(520) are estimated, which are consistent 
with C 2 * 

1. Estimate a ratic L (4 43) /L w (55 0) consistent with 

C 2 by inverting the C 1 algorithm (Sq uat ion ( 1 2 ) with case 
I coefficients) using the C 2 values. 

2. Increase L (670) and decrease L (670) to be consis- 

w a 

tent with the new values of L w (443) . 

3. Recalculate I (443), L (520), and L (550) 

w w w 

4. Recalculate C lr C 2 , K(4 90), K (52 0 ) . 

5. Iterate this procedure until C j and C 2 agree. 

Data acquired aboard the R/V Acania during the Optical 
Dynamics Experiment (ODEX) provide a tentative basis for 
assessing the validity and performance of the above adjust- 
ment algorithm. In Figure 6 values of 1/K(490) calculated 
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from CZCS data, acquired on 16 October 1982, are compared 
with preliminary calculations of l/K(490) from selected C D 3 X 
stations. The transect shewn is along 35N (partition 4) . 
Stations 24 and 25 were occupied 1.5 hours before the Nim- 
bus-7 CZCS observation, and 2 hours after it, respectively. 
Station 21 was occupied S hours, and stations 19 and 13 one 
and two days, respectively, prior to the satellite pass. 
The 1/K(490) values at these stations were calculated from 
the graphical displays of raw irradiance profiles (at a 
wavelength of 490 nm) presented in the preliminary R/V 
ACANIA ODEX CRUISE REPORT (Mueller, Zaneveld and Smith 
1982) . 

Panel 6a compares the CZCS and in situ 1/K (490) values 
before the above adjustment was applied, and figure 6b com- 
pares them after the correction. Agreement in both cases is 
excellent in the transparent waters at stations 21, 24 and 

25: no adjustment for L w (670) was required in this region. 

In the inshore portion of the transect, however, agreement 
is obviously poor before the L (670) adjustment, and much 

W 

improved afterwards. This result is preliminary, and sub- 
ject to possible revision by cognizant ODEX investigators 
when their data have been brought to publishable form. 
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Na verthaiess, the L (670) adjustment algorithm sc overwhelm- 

W 

ingly improves the CZCS estimates of K(490) than its use in 
this thesis project is fully justified and essential. 





DISTANCE (KM) 



Figure 6. Comparison Plots For 1/K (490) Between Track 4 and 
Selected ODEX Stations. Panels (a) and (b) 
Respectively are Before and After Adjustments for 
Non-Zero Values of L (670) . 

w 
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C. DATA STRUCTURE 

Figures 7 through 10 shew the optical depth parameter, 
1/K(490) = Z 90 {490) along each track for the available data 

scenes. (Gordon and HcCluney ( 1975) showed that Z 90 (a) is 
the depth over which 90 percent cf L ( A) is backscattered. ) 
The plots are oriented sc that the ccast is on the right- 
hand side (positive x) , while time of the data scene gees 
from earliest to latest in the positive direction along the 
ordinate cf each figure. The scale of 1/K (490) is in 
meters. 

Chapter II and III give background into the oceanography 
of the region and how that can be related to ocean optical 
parameters. The structures depicted in figures 7 through 10 
will be discussed in terms cf ocean eddy and front visuali- 
zations which result from these relationships. Relatively 
high values of 1/K(490) indicate water with lower concentra- 
tions of chlorophyll and sediment. In general these concen- 
trations may be expected to decrease with distance offshore. 
Abrupt changes in 1/K (490 ) are usually associated with ocean 
frontal structure and eddies. 
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1 . 



Pa:rtitj.on _1 (Zonal Trans sex at 35 5311) 

Th? only data available from 1979 (23 Nov) is from 

winter and shows relatively little structure, (Fig. 7) . 
This image was obtained after the end of the upvellirg sea- 
son and the surface waters were homogeneous to at least 225 
km offshore. 

The 1980 data series shews mere structure. Begin- 
ning on 17 May 80, an eddy of approximately 40 km diameter 
was centered approximately 180 km offshore. Sixteen days 
later the entire track shows several eddy-like features 
ranging in size from 4 km to 20 km. Three days later, on 6 
June 60, the track has lost much of this structure, although 
1/K(490) generally increases in the offshore direction. 
This trend persisted and strengthened slightly through June 
1980. By 1 August 1980, a distinctive pattern had developed 
with nearly uniform turbid waters adjacent to the coast, and 
an abrupt (15 km) frontal transition to much more transpar- 
ent waters at a distance approximately 95 km offshore. This 
pattern is suggestive of the zonal scale of bio-optical 
response to coastal upwelling over a single season. 
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Figure 7. The Optical Depth Parameter, 1/K(490), Across 
Partition 1 (35 53N) 
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ture, except for a weak, clear -water eddy signature 200 km 
offshore. Between 30 September and 16 October, more trans- 
parent (20 to 25 m optical depth) water intruded from off- 
shore to within 190 km of the coast, with a frontal 
transition region of scale approximately 15 km. Over the 
ensuing month, this clear intrusion appeared to evolve into 
a field of less organized, eddy-like anomalies with scales 
ranging from 10 to 50 km. The 2 -dimens ional character of 
this pattern evolution should be studied in a future analy- 
sis of this data set. Such an analysis, which is beyond the 
scope of the present thesis, may indicate whether these 
changes are hesr interpreted as breakdown of a spatially 
continuous intrusion of offshore water, or as simply due to 
advection transporting entirely different water mass 
features into view at this location. 

2. Parti ti on 2 (Zonal Transect at 35 40 N) 

Track two (Fig. 8) is located along latitude 35 
40* N (south of the first track) . 
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Again, the 1979 data set is limited tc one day lane 
in the year and is virtually featureless. 

The 1980 series starts in late spring (6 hay 1990) 

and progresses through late summer (1 August 1980) . The 
same general features found in the northern track are also 
evident in track 2. However, the general appearance of the 
frontal boundary is broken by additional eddy structures. 
In addition, organized ocean frontal structures, which are 
apparently related to similar features in Partition 1 data, 
are here displaced approximately 15 km further offshore. 

The 1982 scenes from partition 2 show the same 
front/eddy development presented in the data from the north- 

* i 

ern transect (partition 1). An offshore frontal boundary 
formed between 30 September 1982 and 16 October 1982, and 
evolved to a less organized pattern of eddy signatures. 

3. Par tit ion 3 (Zonal Transect at 35 22 N) 

Along track 3 (Fig. 9) the features already dis- 
cussed for the previous partitions occurred. An additional 
feature here is a strong clear-water eddy signature approxi- 
mately 220 km offshore, evident on the third of June 1980. 
The data suggest an eddy-like intrusion of transparent off- 
shore water, approximately 40 km in diameter. (A feature of 
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si ailar disaster, but centered approximately 40 km clcssr 
i?. shore appears simultaneously in partition 4 (sea below)}. 

This feature is short-lived , however, as no evidence of it 
remains three days later. This strange phenomenon is 
another candidate for future investigation in a 
2-dimensional analysis. 

In general, frontal boundaries along partition three 
are farther offshore in both the 1980 and 1982 series than 
apparently related frontal signatures in partitions one and 
two. Again, this may suggest the importance of bottom 
steering of the mean flow ever diverging isobaths. 

4. P art ition 4 (Zonal Tran sect at 35 00N) 

The southernmost track at 35 00 N (Fig. 10) has sim- 
ilar patterns, including the offshore displacement of fron- 
tal boundaries. The location of the frontal boundary is at 
its farthest offshore position in this track. 

The transient eddy feature of 3 June 1930, as dis- 
cussed in partition 3, appears here to have a latitudinal 
extent large enough to span partitions 3 and 4 which are 
roughly 45 km apart. This is consistent with the zonal 
scale of the feature (approximately 40 km) . 
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Fiqure 9, The Optical Depth Parameter, 1/K(49Q) , Across 
Partition 3 (3t 22N) 



77 



There also is 3 longitudinal displacement of approximately 
40 km which suggests an eddy of slightly oblong shape that 
roughly parallels the coast. Questions of whether such fea- 
tures are coherent and continuous from one transect to 
another can be easily resolved in future 2 -aimsn sior.al 



anal yses . 

It is reasonable to expect a band of surface water 
with low optical depth to lie adjacent to the coast, due to 
enhanced nearshore biological productivity during the 
upvelling season and due tc northward transport of more tur- 
bid water masses by the Cavidson Inshore Current in winter. 
In contrast, offshore waters tend to be far more transpar- 
ent, at least during the seasons covered by this data set. 
It may be anticipated therefore, that an obvious optical 
front will persistently delineate the boundary between what 
may be classified as nearshore and offshore bio-optical 
regimes, and that intrusions of eddy- like surface water fea- 
tures from one regime to the other will be illuminated by 
optical contrast. Additionally, horizontal gradient struc- 
ture in bio-optical processes will, within each regime, 
often (but not always) accompany the physical structure 
associated with ccean fronts and eddies and produce optical 
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depth structure of similar scales. The scales of previously 
observed structures in the California Current region nay 
thus be expected to be present in the horizontal structure 
of optical depth. It is clear from the above discussion and 
the review of Chapter II that this is in fact the case. In 
the next section these structures will be discussed frcm a 
statistical viewpoint using EOF analysis. 

D. EOF ANALYSIS 

1 • Eigenval ues and D egrees of F re edom 

Table II is a listing cf the eigenvalues for the 
spatial covariance matrix calculated from the data in each 
partition analyzed. The cumulative percentages cf the total 
variance are included. The first ten eigenvalues are listed 
here. In each case, they account for roughly 98 percent of 
the variance. The eigenvalues are presented graphically in 
figures 11 through 14, which include both the eigenvalues 
listed in Table II, and the additional eigenvalues that con- 
tain the nciser-appearina, higher spatial frequencies cf 
variability (and together account for less than 2 percent of 
the total sample variance). We have assumed that 98 percent 
of the variance is an adequate cutoff for calculations . 
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Figure 10. The Optical Depth Parameter, 1/K (490) , Across 
Partition 4 (35 00N) 
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either affected by additional factors 
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2. Data Heccnsn ruction Usi n q Eigenvectors and Pr in cipa l 
ComDcnenos ~ ~ 

Before the eigenvectors and principal components are 
interpreted, how they are combined to reconstruct a particu- 
lar observation is explained. Recall that each eigenvector 
defines a direction of spatial varibility, and that its 
associated principal components represent the amplitude of 
variations in that direction at certain time points. In the 
present context a "direction" takes the form of 1/K(490) 
variations that are coupled at all grid points of the 
domain, and "direction" in this sense may be best 
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TABLE II 





Eigenvalue 


Data for 


Partitions 1 th 


rough 4 




Partiti 


on 1 


Pa r ti 


tion 2 


Order 


Eigenvalue 


Cumulative 


Eigenvalue 


Cumulativ e 




(® 2 ) ~ 


P ercen tage 


2 1 

(m ) 


Percentage 


1 


271. 40 


4 1. 60 


365.90 


43.70 


Ct 


174.90 


68.41 


217.40 


69.67 


3 


72. 17 


79. 47 


99.00 


81 .49 


4 


37.72 


85.25 


4 5 . 60 


86.94 


5 


26. 90 


89. 37 


28.85 


90.39 


6 


17.69 


92.08 


20.39 


92.83 


7 


14.80 


94. 35 


15.12 


94.64 


8 


10.59 


95. 97 


13.63 


96.27 


9 


7. 02 


97. C5 


7.90 


97.21 


10 


5.08 


97. 83 


5.98 


97.92 




Partition 3 


Parti 


tion 4 


Order 


Si gen value 


Cumulative 


Eigenvalue 

2 


Cumulative 




(m ) ♦ 


? ercen tage 


(m ) 


Percentage 


1 


557.40 


54. 18 


273.70 


35.52 


2 


210.60 


74. 65 


183.60 


59.35 


3 


97.76 


84. 15 


11 9.80 


74.90 


4 


58.08 


39 . 79 


44. 14 


80.63 


5 


32.93 
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40.87 


85.93 
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35.96 


90.60 
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13.86 


96. 35 


24.42 


93.77 
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9.48 


97. 27 
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8. 29 
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8.22 


96.58 
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5.74 


98. 64 
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illust 
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r as a 


curve (for 
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transects) or contour plots (for two-dimensional domains). 
Given a temporal mean value for each grid point of a domain, 
the eigenvector multiplied by the principal component for a 
specific time yields a modifier tc the mean signal. 

For example, to view the contributions of the first 
five eigenvectors to the observed signal for 3 June 1980 
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Eigenvalues 



for Partition One 
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EIGENVALUES 




Figure 12. 
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Partition Two 
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Figure 13. 



Eigenvalues 



for Partition Three 
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Figure 14. 
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from partition three of VIS113 
structed curve acres with the 
that a fair reprssen ration can 
the first two eigenvectors. 

3. Mean Structure 



on 25 Jun 1 980 and the reccr.- 
use of five eigenvalues, and 
be reconstructed using onlv 



The mean optical depth (1/k(490) = Z 90 ) transect 

profiles for partitions 1 through 4 are illustrated in Figs. 
16 through 19 and are repeated for ease of comparison in 
Figs. 24 through 27. The mean vector in each track 
represents the tendency of the signal, while the eigenvec- 
tors scaled by the principal components give the perturba- 
tions of the mean. In all four partitions, the mean value 
of optical depth tends tc increase with distance offshore. 
This tendency is expected since the coastal waters should 
contain higher concentrations of sediment and phytoplankton, 
especially during the upwelling season (Traganza, et al. , 
1979). There is a general lack of significant eddy-like 
structure in the mean vectors from all four transects 
(although very lew amplitude perturbations of scale five-to- 
te n km are apparent in the means) . 
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Figure 15. Reconstruction of Optical Depth Transect of 3 

June 1980 (Partition . 3) frca Mean (Solid Curve, 
Panel a) and Successive contri butions of 
Eigenvectors 1 to 5 (Dashed Curves) in Panel a 
to e Respectively. 
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4. Structural Content of Si^en vector s and Princ ipal 
Co mpcrsn ts ” “ ~ ~~ 

The eigenvector discussion involves many intercom- 
parisons of the partitions. Each partition's first ten 
eigenvectors are plotted, (Figs. 16 through 19 and 24 through 
27). To organize the discussion, the first e iaen vector will 
be discussed for all four partitions before proceeding to 
discuss the second, and sc forth. The associated principal 
components are also illustrated, (Figs. 20 through 23 and 28 
th ro ugh 3 1). 

The structure in the first eigenvector of each of 
the partitions is characterized by a band of low variability 
adjacent to the coast, and the structure offshore of that 
band is dominated by a scale extending from there to the 
offshore end of the domain. The "node" marking the onshore 
limit of significant variation in this mode is progressively 
farther offshore, proceeding from the south through the par- 
titions. The "node" of partition 1 begins at approximately 
45 km offshore, and by partition 4, the "node" is 100 km 
offshore. There is a tendency for variance to decrease in 
the amplitude of the first eigenvector as the offshore 
boundary is approached. This may be an artifact of the 
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outer boundary and should be investigated further over a 
larger domain to better estimate the dominant scale. Parti- 
tion 4 , which has the largest spatial extent, shows more and 
larger offshore structure to beyond 180 km. 

The associated principal components, which modify 
the eigenvectors before they are applied to the mean, show 
the time variations. Across the four tracks, the first 

ei genvsct crs/principal components vary in phase with each 
other. In all four transects there is a large difference 
between the first principal component of the only 1979 image 
(early winter) and those from 1980 (early spring) . This 
marked difference is certainly a manifestation of seasonal 
variations in the California Current system (Pavlova, 1966; 
Hickey; 1979) . Most of the first principal components vari- 
ability in all cases is observed in the 1990 series (upwell- 
ing season) , and the record contains relatively little 
variability in the 1982 series (Davidson Current season). 
Coherency of the variations differs from pa rtiti on-to- parti- 
tion with no apparent pattern. The first eigenvector and 
principal component appear to have their foundations in the 
offshore seascnal variation and large scale eddy structure 
that occurs during the upwelling season. In the first mode. 
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the inshore zone influenced by upweiling tends to repair, 
turbid throughout the year, whereas the dominant variations 
in optical depth occur offshore of the upweiling zone. 

The shapes of second eigenvectors from the four par- 
titions are similar, but vary in an oscillatory fashion from 
the northern partition to the southern partition. The first 
partition (northern) shows negative values beyond approxi- 
mately 180 kra offshore, and then smali amplitude positive 
values from there to the coast . Partition two and three 
depict a mirror image pattern to that of the first parti- 
tion. Partition four shews much the same pattern as parti- 
tion one. The phase relation in the principal components 
shows no pattern between partitions one and two, but the 
series for partitions three and four both suggest phase 
reversal from the first partition. This negati ve-to-posi- 
ti ve-back- to-negati ve pattern of behavior weakly suggests a 
wave-like meridional oscillatory structure, with an offshore 
peak (180 km offshore) in the vicinity of partitions two and 
three. Resolution of this meridional characteristic feature 
will require a 2-dimensional analysis. The distance of sep- 
aration of the partitions suggests a wavelength of the order 
of 120 km. Again, the majority of the variability occurs 
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Figure 16. Hear and Eigenvectors 1 to 5 for Partition One. 
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17. Mean and Eigenvectors 1 to 5 for Partition Two. 
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Figure 18. Mean and Eigenvectors 1 to 5 for Partition 
Three. 
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Mean and Eigenvectors 1 to 5 for Partition Four. 
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Figure 20. Principal Components 1 to 5 for Partition One. 
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Figure 21. Principal Components 1 to 5 for Partition Two. 
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Four. 



during the I960 series. More nearshore structures are 
apparent in this eigenvector, as is an increase in nearshore 
variability, as compared to the first eigenvector. 

The third eigenvector has a similar behavior for the 
first two partitions. The perturbations are of roughly the 
same spatial scale (45 km) and appear to be in phase. How- 
ever, the third and fourth partitions show an opposite 
behavior in the far offshore region (beyond 180 km). Numer- 
ous smaller scale features (of the order of 10 km or less) 
are apparent in this eigenvector. In general, higher spa- 
tial frequencies become increasingly important in higher 
order eigenvectors. The principal components show an 
increase in the variability of the 1982 series with the wide 
range of variability still present in the 1980 series. This 
eigenvector shows the largest nearshore amplitudes of all 
the eigenvectors, which suggests it may be closely linked to 
the nearshore structure cf upwelling. There is little sug- 
gestion of a temporal relation evident in the third princi- 
pal components of the four tracks. 

The fourth eigenvector shows an increase in fre- 
quency (decrease in wavelength) cf the represented variabil- 
ity scales. Features range in size from 18 - 45 km, with 
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numerous smaller scale perturbations. The average 

wavelength and range of variability is approximately the 
same for all four partitions. The principal components for 
pa rtitions one , three and four qualitatively suggest coher- 
ency. However, eigenvector behavicr is opposite in parti- 
tion three and similar in partition four when compared to 
partition one. A prominent feature approximately 180 km 
offshore in the fourth eigenvector of partition one, seems 
to be shifted outward to 210 km offshore in partition four. 
This time/space relationship between the structures of par- 
tition one and four again suggests a meridional oscillation 
worthy of future investigation through a 2-dimensional 
analysis. 

The fifth eigenvector shows features of scale that 
range from 5-42 km. The much smaller features (less than 5 
km) are not dealt with as they are essentially part of the 
background noise expected in any natural system. Little can 
be said about the correspondence of the four partitions with 
just a visual inspection. However, nearshore structure 
appears in the eigenvectors with more variability than in 
previous eigenvectors of the 1982 series. The 1980 series 
still demonstrates the largest overall variability, and the 
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large inter annual differences between the one scene in 1979 
and the 1980 scenes is still evident. The principle compo- 
nents of partition one and four show excellent agreement in 
amplitude and phase for the first five time points, but as 
structural variability decreases with time, so does corre- 
spondence. The quantitative correspondence of these varia- 
tions is beyond the scope of this thesis, but it should be 
investigated in future analyses of this data set. 

The sixth through the tenth eigenvectors are charac- 
terized by variation of such high frequency that little can 
be said of the relationships between the partitions. Scales 
of structural features in these eigenvectors range from 1 to 
35 km, with nc suggestion of a temporal relation between 

I 

partitions. The principal components show that the 1980 
data again dominates the variance, but at these higher fre- 
quencies the increased contribution cf the 1982 data to the 
total variance is very apparent. Eecause of this disorgan- 
ized structure, detailed interpretation is not attempted for 
eigenvectors cf order greater than 5. 
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Figure 24. Mean and Eigenvectors 6 tc 13 for Partition One. 
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Figure 25. Mean and Eigenvectors 6 tc 13 for Partition Two. 
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gure 26. Mean and Eigenvectors 6 tc 10 for Partition 
Three. 
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Figure 27. Mean and Eigenvectors 6 tc 10 for Partition 
Four. 
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. Principal Components 6 to 10 for Partition Two. 
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gure 30. Principal Components 6 to 10 for Partition 
Three. 
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5. The Join ing Of Tvc Partit ions 

Only a cursory treatment of the joining function 
results is presented. Partition one and two were examined 
and analyzed as per the joining function development given 
in Chapter IV. The analysis was based on performing the EOF 
analysis on partitions one and two combined, and then com- 
paring this result to these obtained from the separata EOF 
analyses that were joined. It was found that the joining 
function principle components were within .0001 m of these 
computed using the two partitions as one data set. Further- 
more, using ten degrees of freedom, it was found that 

T " 6 

^j 1 J = I (within 2.5 X 10 ) (51) 

and 

T • . ' 6 

E 1 E = I (within 8.0 X 10 ) (52) • 

'Xt Oj Oj 



This result is based on the orthogonality of the eigenvec- 
tors as J is the matrix of eigenvectors of the covariance 
matrix of principal components for the two partitions com- 
bined. 

Additionally, a comparison of the principal compo- 
nents yielded (For a joint sample size of 17) 




- 3 

(within 1.4X10 ) (53) 
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L 

% 



Representing the principal components of the joining func- 
tions 



with Y , it follows that 



(l6 ) 



Y Y = L 

'x, a- a. 



(within 5.0X10 ) ( 54) 



Again the small difference demonstrates the utility of the 
joining process. 

The eigenvalues obtained by thus joining the eigen- 
vectors and principal components from partitions 1 and 2 are 
given in Table III, together with fractions of total sample 
variance. At face value, the first ten eigenvalues account 
for 98.55% of total sample variance. Rscall however, that 
the input data were represented in truncated form, using 
only the first ten principal components from each of the 
partitions. This original a p pro ximati on retained only 
97.89% of the total variance computed from the original 
data, hence, it is necessary to adjust the apparent trunca- 
tion cf the joined result accordingly. The results of this 
adjustment are given in the third column of Table III and 
show that assuming that only the first ten eigenvectors are 
significant actually leads to a truncation to 96.50% of the 
total sample variance. While 3.5% precision is an accepta- 
ble level cf approximation for most problems in geophysical 
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data interpretation 



the effects of 



successive 



must fce giver, careful attention when applying the parti- 
tioned method to EOF analysis. 

This example indicates strongly that the partition- 
ing approach to EOF analysis can provide computationally 
acceptable results when applied tc satellite image data. 
The example also emphasizes that proper care must be given 
to controlling successive truncation in the partition join- 
ing process. It is left to future projects to investigate 
questions such as joining partitions on the basis of par- 
tially intersecting samples to provide optimal functions 
for interpolating satellite data into cloudy regions, and to 
interpretation of partition joining functions to illuminate 
spatial correlations between locally important structures 
(e.g., topographically generated mesoscale eddies) and the 
dominant structure of the overall domain (e.g., that associ- 
ated with -he evolution of the synoptic scale upwelling 
front over the continental slope and shelf over the course 
of the upwelling season) . It is questions of this kind that 
address the ultimate utility of partitioned EOF analysis. 
The present effort is limited tc preliminary work to estab- 
lish foundations of feasibility and procedural constraints. 
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TABLE III 



Eigenvalue Data for Joining Process 



gen val ue 


Percenta ge 


Cumulative 


(m ) 


of Variance 


Pa rcan tage 


594 .00 


40.74 


40.74 


401.60 


27.54 


68. 28 


17 1.70 


11.78 


80.06 


67.35 


4.62 


94.68 


6 3.86 


4.38 


89.06 


47.57 


3.26 


92.32 


34.34 


2.36 


94.68 


23.17 


1.59 


96.27 


20.52 


1.41 


97.68 


12.70 


0.87 


98.55 
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VI . 



2£SCiiSSI0N AND conclusions 

Zonal transects of optical depth (1/k(490) m) measured 
with the Nimbus-7 CZC3 have been analysed to investigate 
bio-optical structure over the continental shelf and slope 
off central California. Samples of cloud free data were 
selected and processed for latitudes 35-53N, 35-4QN, 35-22N 

and 35-00N. The data were observed in 1979, 1980, and 1982 
during the months Hay through November. The zonal structure 
in these samples was analysed using EOF's computed sepa- 
rately for each section. Meridional variance structure was 
analysed only qualitatively through inspection of similari- 
ties in features contained in EOF's of the different tran- 
sects and in the temporal sequences of associated principal 
components. Finally, the computational feasibility of 
applying partitioned EOF analysis methods to this type cf 
data was investigated by joining the EOF's of the two north- 
ernmost transects to form estimates of the EOF's of the com- 
bined spatial domain. 

The first eigenvectors for four zonal transects of opti- 



cal depth 1/(kU90) each contained dominant scales of order 



290 km or greater, and accounted for between 35 and 54 per- 
cent cf the total variance. They are each also character- 
ized fcv a band cf low variability in optical depth in the 
inshore region influenced by upweiling and the Davidson 
Inshcre Current. This band is confined within 45 km of the 
coast at 35 53’ N, and mcnoton ica lly broadens to approxi- 

mately 100 km at 35N latitude. This behavior is possibly 
related to the broadening of the continental shelf and slope 
with longshore distance scuth of hcnterey. Hurlburt (1979) 
showed that the topographic beta effect plays a fundamental 
role in the dynamics associated with mesoscale (order ICO 
km) longshore variations in topography by affecting the 
strength cf the longshore flow. Also, the influence of 
topography can produce barctropic flew beyond its immediate 
vicinity. For mesoscale variations in coastline geometry, 
the coastal currents and the patterns of vertical motion 
tend to follow the coastline, but net with uniform strength. 
Coastal current widths tend to be narrower than the scale of 
coastline variability. In these terms, the meridional vari- 
ation in scales present in the first EOF' s are qualitatively 
consistent with the longshore variations in bathymetry of 
the study domain. 
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The second eigenvectors account for zonal structure with 
dominant scale cf order 120 km, and with nearly uniform 
amplitude from the ccast tc a node approximately 150 km off- 
shore in partitions 1, 2, and 4. The second eigenvector for 
partition 3 (35 20N) is anomalous in that it is dominated by 
a zonal waveform with nodes spaced at approximately 80 km, 
or roughly half the dominant scale of its counterparts. The 
reason for this behavior should be investigated. 

The third eigenvectors are dominated by scales ranging 
from approximately 60 to 100 km (between nodes) . The shapes 
and scales vary more strongly from partition-to-part ition 
than was the case with the first two eigenvectors. 

Across each transect, zonal features with wavelengths 
100 km and greater appear. The suggestion of an oscillatory 
behavior in the meridional direction needs to be studied 
further. Resolving such a feature requires a more detailed 
study involving a 2- di mansion al analysis cf the study 
domain. 

The large eddy field associated with the shoreward 
boundary of the California Current was observed in the data 
sat. The scales of this eddy field were of the same 
magnitude as the spatial scales employed in the partitions. 
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this 



This necessitated placing the partition boundary within 
eddy field arc cutting away seme of the features. The 
sor.al development of a synoptic scale upwelling front off 
the California coast is strongly suggested in the data and 
its eigenvectors. The smaller eddies associated with this 
pattern ranged from 5 to 100 1cm in scale. 

The convergence of the eigenvalues to roughly 98 percent 
of the variance after the tenth value was of particular 
interest. This was true for all four partitions and 
although this is not an overwhelming reduction in the 
degrees of freedom of the initial system, it is significant. 

Satellite images, and ether fields of oceanic and atmos- 
pheric variables, provide massive data sets. Large amounts 
of computer time must often be expended for processing these 
data sets at even relatively primitive levels. Analyses and 
interpretations are, morever, made difficult by the sheer 
volume of data. EOF analysis provides a viable method for 
mathematically representing satellite data fields in a com- 
pact and easily manipulated form. Data transformed using 
EOE’s illuminates, and facilitates analysis of, the time and 
space scales associated with a given variable over the 
domain; the present study has exercised this attribute of 
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EOF's on a descriptive level. In addition, the contact 
principal component representation of satellite images pro- 
vides an efficient form for analysing the response of spa- 
tial structure in, for example, optical depth to forcing by 
wind stress and currents, acting through a bio-optical 
model; this is a logical avenue fcr future research to build 
on the present results. 

Considering purely computational aspects of EOFs, the 
well-known symmetry of eigenvector solutions in the time and 
space domain can be used tc great advantage in the analysis 
of satellite image data. The number of spatial grid points 
in even the single trackline partitions of the present 
study yield large, but computationally tractable, scatter 
matrices. The larger arrays associated with 2-dimensicnal 
area partitions, each with several hundred grid points, will 
clearly exceed sizes admitting direct computation of spatial 
EOFs. The linear algebra and scalings involved in using the 
smaller time domain scatter matrix for computation of space 
domain EOF's is reviewed in Appendix C. 

The partitioned method of EOF analysis illuminates cor- 
relations between variability in spatially separate sub-re- 
gions. The present results demonstrate the computational 
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faasiblity cf this piecewise approach when applied to CZCS 
optical depth data. There is every reason to believe -hat 
the method may be equally well applied to other CZCS parame- 
ters and tc infrared imagery of SSI. Further work in this 
area should aim to determine whether the joining functions 
linking EOF's from separate domains are sufficiently sta- 
tionary tc provide a basis for optimally interpolating sat- 
ellite image data of these types over cloud-covered areas of 
a particular day's image. Other applications to be explored 
include determination of the extent to which correlations 
between 3 -dimensional in situ data and 2-dimensional satel- 
lite data in small sub-regions may be extended to other 
parts of the larger domain covered by satellite data alone. 
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APPEND I X A 



SATELLITE DMA PROCESSING METHODS 



A. INTRODUCTION 

Data processing was divided into three major levels. 
Lavel-I processing includes all steps required to take the 
original data tape to a Level-I tape. Level-II processing 
includes all steps between a Level-I tape and a Level-II 
tape. Level-Ill processing includes the steps involved to 
take the Level-II output to a usable form. The following 
sections briefly describe the steps involved in the three 
levels of processing. 

Computer hardware utilized was that resident at the 
Naval Postgraduate School, Monterey, California. The main 
frame computer used was the I3M 3033A? while the mini-com- 
puter used was the AppleMI. Computer software referred to 
in this section is either a system utility resident to the 
IBM system or a locally generated program. Documentation of 
the locally generated programs can be obtained from: 

Dr. J. L. Mueller (Cede 68My) 

Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93943 
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The software involved has many system-depet 
well as features inserted for convisnence. 
programs cn other systems are cautioned to 
mentation carefully prior to attempting 
software . 



Users of these 
review the dccu- 
to transfer the 



B. LEVEL-I FfiOCESSI NG 

figure 32 is a schematic diagram illustrating the 
processing steps for Level-I and should be referred to 
throughout this discussion. The master tape (raw satellite 
data) was obtained from the Scripps Institution of Oceanog- 
raphy, San Diego, California. Table 17 gives a summary of 
the master tapes utilized in this study. The data were in 
the form of a ^standard magnetic tape in a binary format with 
6250 tits per inch (BPI) . The tapes wera originally created 
using a Hewlett Packard (HE) - 3000 which has a characteris- 
tic high crdar, low order bit arrangement opposite to the 
IBM system. Therefore, before using this raw data in the 

13 M 3033AP, it had to undergo a byte swap routine. This 
byte swap was accomplished when the unformatted working bac- 
kup tape was made using local pregram VISBKV. After the 
unformatted backup tape is made a variable blocked spanned 
(VBS) format tape is produced using the system utility 
IEBGENER. 
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Figure 32. 



Level-I Processing Schematic Diagram 
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TABLE IV 



'Tape 

Da signs ticn 
VISO 17 
VI S 0 3 2 
VI S040 
VIS094 
VI SO 95 
VI S097 
VIS 1 04 
VI SI 05 
VIS1 06 
VI S 1 17 
VI S 1 1 1 
VIS1 12 
VI SI 13 
VIS1 26 
AR0000 
AR 26 42 
AR2668 
AR 2 6 85 
AR 26 86 
AR 26 9 1 
AR 26 93 
AR 27 04 Nimbus 



its Data Tape 
Da-s 

(CZCS) 16 OCT 

(CZCS) 12 NOV 

(CZCS) 23 NOV 

(CZCS) 6 MAY 

(CZCS) 5 MAY 

(CZCS) 17 MAY 

(CZCS) 3 JUN 

(CZCS) 6 JUN 

(CZCS) 7 JUN 

(CZCS) 12 JUN 

(CZCS) 23 JUN 

(CZCS) 24 JUN 

(CZCS) 25 JUN 

(CZCS) 1 AUG 

(CZCS) 30 SEP 

(CZCS) 5 OCT 

(CZCS) 16 OCT 

(CZCS) 27 OCT 

(CZCS) 28 OCT 

(CZCS) 1 NOV 

(CZCS) 3 NOV 

7 (CZCS) 14 NOV 





S3 T avail 


1979 


no 


1 979 


nc 


1979 


no 


1980 


yes 


1980 


yes 


1980 


yes 


1980 


yes 


1980 


yes 


1980 


yes 


1 980 


yes 


1980 


yes 


1 980 


yes 


1 980 


yes 


1980 


yes 


1 982 


no 


1982 


no 


1982 


no 


1 982 


no 


1982 


no 


1982 


no 


1982 


no 


1932 


no 



Sat ail 

So u r ce 

Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
Nimbus 7 
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This format is used in conjunction with unformatted read 
statements to minimize computer time. These two cooled 



of ~ h 3 



tapes serve as the working tapes for the remainder 
Level-I processing, and the master tape is archived. 

Using the VBS formatted tape, a Versatac plotter grays- 
cale is produced using local program ZIPSIO. This program 
also unpacks the 2vent file (hereafter referred to as the 
2-file) and writes it tc a storage disk. The grayscale 

depicts the satellite pass in picture form for hand analysis 
of landmarks. line numbers and pixel numbers are taken off 
the grayscale for clear, cloud-free landmarks. These values 
are entered into local program ZIPPIC to generate a 
' P ICPRINT ' . This is a matrix of radiance value centered on 
the individual landmarks line and pixel number. These 
PICPRINTS are then contoured by hand (using a threshold 
value of 18 counts for land or clouds) tc determine an exact 
time and pixel number fcr the landmark. The landmark's 
latitude, longitude, line number ar.d pixel number with addi- 
tional housekeeping data are entered into local program 
CZCSNAV on the Apple II. This program is interactive and 
prompts for necessary inputs. Additionally this program 
adjusts roll, pitch and yaw to reduce the root mean squared 
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distance error in the navigation problem. The mean ms 
value obtained for all the adjusted, utilized data was 
approximately 1.09 n.mi.. The final product of this step 
generates a set of navigation parameters that are used to 
generate a navigation matrix. This step is accomplished 
using local program CZCSUAV2 to generate the navigation 
matrix and NAVDUMP to write the navigation matrix (here- 
after referred to as the G-file) to a temporary formatted 
tape. 



The E-file is copied from its temporary disk storage to 
the temporary tape as the G-file. Additionally the Data 
file (hereafter referred to as the L-file) is first reversed 
from its fccttcm-to- top orientation to a top- to- bottom orien- 
tation using local program BACKWARD. This program also puts 
the L-file tc the previously mentioned temporary storage 
tape. Finally, these files on the temporary tape are copied 
to a Level-I tape using the system utility IEBGENER. The 
only difference between the temporary tape and the final 
Level-I tape is that the L-file is copied into an unformat- 
ted file which will aid in the speed of further processing. 
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C. LEVEL -II r 50 CES S I NG 

Figure 33 is a schematic diagram illustrating the 
processing steps for Level-II. The Level-I taps generated 
by the steps discussed in the previous section is the input 
tape for this processing. Only the L-file is affected by 
the Level-II processing as the E-file and G-file are copied 
straight to the Level-II tape using the system utility 
I53GENER. The L-file is used to generate output for calcu- 
lating the proper values cf the Angstrom coefficient for 
each scene. This is done using local program CZ PAR M S 2 and 
an assumed value for the Angstrom coefficient. Chapter III 
Section C.2. discusses the importance and method of finding 
these values. Next, the computed Angstrom coefficients with 
the L-file are rerun through CZPARNS2 to regenerate the 
L-file. This regeneration involves talcing the raw counts of 
each channel and applying the bio-optic algorithms discussed 
in Chapter III tc produce values fcr chlorophyll and K. At 
this point the adjustment algorithm discussed in Chapter V 
has not been applied. Analyses cf the initial Level-II out- 
put precipitated the need for the corrective algorithm, 
which was then applied during Level-Ill processing. 
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33. Level-II Processing Schematic Diagram 
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D. L ZV EL- III PROCESSING 

Figure 34 illustrates the Level-Ill processing steps. 
The Level-III processing basically takes the data obtained 
in Level-II and marries it no the navigation matrix gener- 
ated during Level-I processing. Using the four designated 
tracks (Fig. 5 ), the coastal starting points from each 
track were entered into local program TED DUMP to provide the 
navigation block of the track origin. The G-file contains 
data (in latitude and longitude values) every sixteenth 
pixel and sixteenth line. Once the origin block is estab- 
lished the exact line and pixel was interpolated using local 
program FINDPIX on the Apple II. With this starting point 
local program DATA4 was entered to generate every 1 km along 
each track an associated line and pixel number which was 
then converted into the appropriate data values. This out- 
put was written to storage for later processing. It was 
here that the adjustment algorithm was applied, producing 
the final version of the data in a navigated form. 
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APPENDIX 3 



DATA CONDITIONING 

To apply the Level-Ill data tc the analytic techniques 
certain conditioning steps were necessary prior to begin- 
ning. Much of the conditioning applied to the data was 
dependent on the data itself as to its completeness and 
behavior. This discussion focusses on the steps necessary 
prior to using the EOF analysis techniques. 

Figure 35 depicts the steps involved in this discussion 
and should be referred to as a guide. First the Level-Ill 
data for the four tracks and twenty-two scenes were 
extracted using local program PACKJOB. Twenty-two files 
each contained the data for the four tracks for each partic- 
ular scene. These data were plotted using local program 
PARPLOT and the DISSPLA utilities resident on the 
IBM-3033AP. The format of the plot was chosen to give an 
indication of either good data or bad data with no struc- 
ture. This plot was used to decide on the partitioning 
scheme. Four partit ions_were selected and their details are 
listed in Table V. 
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TABLE V 

Partition dimensions 



Track 


Part it 


ion 


Min Grid 


Max Grid 


N t i a e 


Ns pace 


No . 


Mo. 












1 


1 




1 £0 


410 


17 


231 


2 


2 
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17 


26 1 


3 


3 




200 


467 


16 


268 


4 


4 




2 CO 


485 


17 


2 86 


Figur 


es 36 


thro 


ugh 39 show 


the plots 


generat 


sd by p 


plot and 


the par tit 


ioning given 


in Table 


IV. The 


ob j ect 



of the partitioning was tc find the most complete data over 
time and space possible given ten samples. 

Eata from these four partitions was then entered into 
local program PARDAT. This program applied most of the con- 
ditioning to the data set. Only the K(490) data was uti- 
lized from this point on although this program could be 
easily altered to focus on another optical parameter. The 

data were searched to find good points and bad points and a 
control arrangement for later use was made. The raw data 
were scaled and inverted to produce 1/K(490) values in 
meters. Next the data were averaged by every fourth point 
to smooth out noise features. At this point data strings 
with gaps existed for each applicable scene. Next a linear 
interpolating routine was applied to obtain continuous data 
at each time point. Finally, the dara for all scenes and 
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tracks were combined and written to disk in a single data 
file. To this pcint the conditioning applied has consisted 
of partitioning the data into four partitions, rejecting 
incomplete scenes, scaling the K (U90) values, averaging the 
data by every four values, and applying linear interpolation 
to fill in the remaining gaps. 

The conditioned data was then plotted using local pro- 
gram PARPLT and the DISSELA system utilities. The plots 
generated are figures 7 through 10 and were used in Chapter 
V section B to discuss the data and its relationship to the 
regional oceanography. 

The final steps of the data conditioning involved appli- 
cations of the EOF analysis techniques. Local program PEOF1 
produced eigenvalues, eigenvectors, and principle components 
for each partition and plotted the output. Figures 11 
through 14 and 16 through 31 are the plots produced. The 
eigenvalues, eigenvectors and principle components were all 
written to disk for later use. The local program JEOF2 was 
designed to produce the joining function that related parti- 
tion one to partition two. 
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Figure 36. Partitioning Scheme for Track One (35 53 N) 
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Figure 37. Partitioning Scheme for Track Two (35 40 N) 
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Figure 38. Partitioning Scheme for Track Three (35 22 N) 
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Figure 39. Partitioning Scheme fcr Track. Four (35 53 N) 
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APPENDIX C 



EOF PROCESSING 

Ths desired EOF • s are those in the space domain, which 
for satellite data, is dimensioned much larger than the time 
dimension. It is possible tc significantly expedite compu- 
tations by computing eigenvalues, eigenvectors and principal 
components using the smaller covariance matrix of the time 
domain, and to then scale these results to obtain the eigen- 
vectors and principal components in the space domain. The 
algebraic basis for this approach is reviewed in this 
appendix . 

Consider the following convention for dimension nota- 
tions 

Space . . . . n =1,...,N 

Time . . . . m = 1 , . . . , a 

EOF (order) . . . . k = 1,...,K 

where K< min(a,N). As before the raw data matrix is given 
by 

F» = f f * 1 (55) 

% L ran J 

of M rows by N columns. The sample mean in the space domain 
is given by 
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M 



f = M 



E 



i.i -- 1 



mn 



{ 56 } 



of dimension 1 X N. The centered data matrix in the space 
domain is given by 



F = 



f ' - f 

- n 



(57) 



From Chapter IV the sample space covariance matrix (of size 
N X N) is 



1 



S = M-l F F 



(58) 



and the sample time covariance matrix (centered in space and 
of size M X M) is now defined as 



1 



T = M-l F F T 



(59) 



If X, and e, are an associated eigenvalue and eigenvector 
k -k 



respectively cf S, then 

v/> 



S e . - X . e 7 ^ 0 

'v -k k -k 



(60) 



si nee 



X^ ^ 0 and f 0 



(61) 



Then 



F e, t ° 
'v -k 



(62) 



and therefore 
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—k- F F T F e. = X, F e, 
M_ l ^ >\j - k k ^ - iC 



(635 



(noting the commutative property of X^ and £ ). Thus,X k is 
an eigenvalue of T with associated eigenvector (in the time 
domain) 



a, = F e 
-k 'u -k 



<SU) 



where a , is the k'th principal component (in the space 

ink 

domain) at time point m. 

We now consider the eicenvalues x, and eigenvectors 

k -k 



of T, normalized such that 

n . 7 






M 



E z . 5 . = 5 . . 

mi mj lj 



(65) 



m = 1 

with time domain principal components 



5k * E T - 5 k 



k = 1 , . . . , K 



( 66 ) 



where the ' signif icance order K is selected either to retain 
an arbitrary traction of the total sample variance, or 
according to cr.e of the more objective statistical selection 
rules discussed by Preisendor fer , et al. (1981). 

The principal components in the time and space domain 
have the property that 




(M- 1 ) X. 



5 . . 

ij 



( 67 ) 
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Recall that F was centered on f an the space domain, and 
chat both S ar.d T were therefore normalized by the factor 

(1/M-1) • 

Ey analogy tc equation (63) 



ITT T 

— F F F 5, = >•- F 5, = 

M- 1 ^ ^ ^ - k tc a, - k k-k 



(63) 



where 5 , is now seen to be an eigenvactor of s, but of 

” rC 

leng th 



^k 



= / (M - 1) X. 



(69) 



rather than being an orthcnormal eigenvector of length 1 



The two are thus related by 

c. 



-k /(M - 1 ) a 



(70) 



(principal components in the time domain scale to orthonor- 
mal eigenvectors in the space domain), 



and by similar arguments 

- ('(m - m k ) c k < 71 > 

(orthonormal eigenvectors in the time domain scale to prin- 
cipal components in the space domain) . 
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